"0","inds <- data.frame(seq(1,100))
colnames(inds) <- c(\"N\")

EsN.calc <- function(df) {
  rarefied <- rarefy(sampled.inds,df$N, se = TRUE, MARGIN = 1)
  EsN <- rarefied[1]
  Es.se <- rarefied[2]
  outs <- data.frame(EsN,Es.se)
  return(outs)
}

Rarefied.richness <- ddply(inds,.(N),each(EsN.calc))"
